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Abstract  —  This  paper  presents  a  linear  quadratic  regulator  approach  to  the  robust  stabilisation, 
robust  performance,  and  disturbance  attenuation  of  uncertain  linear  systems.  The  state-feedback 
designed  systems  provide  both  robust  stability  with  optimal  performance  and  disturbance  attenuation 
with  H oo-norm  bounds.  The  proposed  approach  can  be  applied  to  matched  and/or  mismatched 
uncertain  linear  systems.  For  a  matched  uncertain  linear  system,  it  is  shown  that  the  disturbance- 
attenuation  robust-stabilizing  controllers  with  or  without  optimal  performance  always  exist  and  can 
be  easily  determined  without  searching;  whereas,  for  a  mismatched  uncertain  linear  system,  the 
introduced  tuning  parameters  greatly  enhance  the  flexibility  of  finding  the  disturbance-attenuation 
robusUstabilizing  controllers. 


•  • 


1.  INTRODUCTION 


The  problems  of  robust  stabilization,  robust  performance,  and  disturbance  attenuation  of  un¬ 
certain  linear  systems  have  drawn  much  attention  recently.  Nonlinear  robust  control  laws  that 
stabilize  uncertain  linear  systems  satisfying  matching  conditions  were  developed  by  Leitmann  [1]. 
Feedback  control  designs  based  on  the  Algebraic  Riccati  Equation  (ARE),  which  adjust  a  scalar 
to  achieve  stabilization  of  the  systems  with  uncertainty  parameters  bounded  by  constraint  sets, 
were  derived  by  Petersen  and  Hollot  [2],  Petersen  [3],  Schmitendorf  [4],  and  Khargonekar  et  al.  [5]. 
These  approaches  have  generally  utilized  the  concept  that  a  given  ARE-based  control  law  guar¬ 
antees  the  existence  of  a  quadratic  Lyapunov  function  (and  hence,  stability)  for  the  closed-loop 
uncertain  linear  system.  Also,  other  recent  research  attention,  e.g.,  Bernstein  and  Haddad  [6], 
Doyle  et  al.  [7],  Glover  and  Doyle  [8],  and  Petersen  [9],  has  been  given  to  the  ARE-based  control 
designs  which  stabilize  a  nominal  system  and  reduce  the  effect  of  disturbances  on  the  output  to 
a  prespecified  level.  More  recently,  Veillette  et  al.  [10]  has  proposed  an  ARE-based  design  which 
not  only  robustly  stabilizes  an  uncertain  linear  system  with  the  structured  uncertainty  in  the 
system  matrix,  but  also  provides  disturbance  attenuation  with  a  robust  H^-notm  bound. 

In  this  paper,  based  on  linear  quadratic  regulator  theory  and  Lyapunov  stability  theory, 
we  develop  linear  state-feedback  control  laws  for  robust  stabilization,  robust  performance,  and 
disturbance  attenuation  of  a  given  uncertain  linear  system  with  the  uncertainties  existing  both 
in  the  system  matrix  and  the  input  matrix.  The  proposed  design  procedures  can  be  applied  to 
both  matched  and  mismatched  systems.  The  paper  is  organized  as  follows.  First,  the  matching 
conditions  for  uncertain  linear  systems  to  be  stabilized  with  prespecified  disturbance  attenuation 
level  are  defined  in  Section  2.  It  is  shown  that  many  dynamic  systems,  described  by  second- 
order  monic  vector  differential  equations,  often  satisfy  these  matching  conditions.  Next,  linear 

•This  work  w«s  supported  in  part  by  the  U.S.  Army  Research  Office,  under  Contract  DAAL-03-91-G0106, 
and  NASA-Johnson  Space  Center,  under  GRANT  NAG  9-380, 
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robust-stabilizing  controllers  which  provide  disturbance  attenuation  and  optimal  performance 
for  matched  systems  with  norm-bounded  or  structured  uncertainty  matrices  are  developed  in 
Section  3.  Also,  it  is  shown  that  linear  disturbance-attenuation  robust-stabilizing  controllers  with 
optimal  performance  for  matched  systems  always  exist  and  can  be  easily  determined  without 
searching.  Then,  in  order  to  achieve  the  robust  stabilization  and  disturbance  attenuation  of 
mismatched  systems  with  norm-bounded  or  structured  uncertainty  matrices,  alternative  linear 
distrurbance-attenuation  robust-stabilizing  controllers  are  proposed  in  Section  4.  To  demonstrate 
the  proposed  methods,  two  examples  are  illustrated  in  Section  5,  and  the  results  are  summarized 
in  the  conclusion  in  Section  6. 

2.  NOMENCLATURE,  SYSTEMS,  AND  DEFINITIONS 


Nomenclature 


^nui(Af) 

l|M|| 

/ 

0 

Af  >(>)  0 
M  < (<) 0 
P>(>)« 
P<(<)Q 


maximum  singular  value  of  a  matrix  M\ 
minimum  singular  value  of  a  matrix  M; 

matrix  norm,  j|M||  =  Om»i(Af)  = 
identity  matrix  of  appropriate  dimension; 
null  matrix  of  appropriate  dimension; 
matrix  M  is  symmetric  positive  (semi)definite; 
matrix  M  is  symmetric  negative  (semi)definite; 
means  P  —  <9  >  (>)  0; 
means  P  —  Q  <  (<)  0. 


Consider  the  uncertain  linear  system 


i(t)  =  (A  +  AA)x(f)  +  (B  +  AB)u(t)  +  Dw(t), 
y(t)  =  C  x(t), 


(la) 

(lb) 


where  x(<)  g  Kn  is  the  state,  u(f)  g  Hm  is  the  control,  u>(<)  g  H1  is  the  disturbance,  y(t)  g  TV  is 
the  controlled  output,  A  g  Hnxn,  B  g  Knxm,  D  g  Rnxi,  and  C  g  Wxn  are  the  nominal  system 
matrix,  input  matrix,  disturbance  matrix,  and  output  matrix,  respectively,  and  A  A  and  A  B  are 
the  associated  uncertainty  matrices  of  appropriate  dimensions  with  respect  to  A  and  B.  Note  that 
the  uncertainty  matrices  AA  and  A B  can  be  time-varying.  We  assume  that  the  nominal  system 
(A,  B)  is  controllable.  Our  objective  is  to  design  a  linear  state-feedback  control  law  u(t)  =  K  x(t) 
such  that  the  resulting  closed-loop  system  matrix  .4C  =  A  +  A  A  +  (B  +  A  B)  K  is  asymptotically 
stable,  and  the  resulting  closed-loop  system  is  optimal  with  respect  to  a  performance  index, 
and  the  Hoo-norm  of  the  closed-loop  transfer  functon  matrix  H(s)  =  C(sl  —  Ae)_1  D  from  the 
disturbance  input  w(t)  to  the  output  y(t)  is  less  than  or  equal  to  some  prespecified  disturbance- 
attenuation  value  6,  i.e.,  f/T(-jw)  A/(jw)  <  62  I  for  all  w  g  H. 

To  proceed  with  the  derivation  for  such  a  control  law,  we  need  to  consider  two  classes  of 
uncertain  linear  systems  which  are  matched  and  mismatched.  The  system  in  (1)  is  called  a 
matched  uncertain  linear  system  if  there  exist  matrices  E  g  7?mxn,  F  g  7 Zmxm,  and  G  g  Hmxi 
such  that 


(i)  AA  =  BE, 

(ii)  AB  —  B  F,  and  ||F||  <  1  or  2/  +  F  +  FT  >  0,  and 

(iii)  D  =  BG. 

The  matching  conditions  (i)  and  (ii)  constitute  sufficient  conditions  [1]  for  the  system  to  be 
stabilizable.  We  shall  show  that  the  uncertain  linear  system  is,  in  fact,  linearly  stabilizable  with 
any  disturbance  attenuation  6  >  0  if  it  satisfies  conditions  (i— iii). 

It  is  important  to  note  that  a  dynamical  system  [11]  which  can  be  modeled  by  a  second-order 
monic  vector  differential  equation  is  often  a  matched  system.  This  fact  can  be  verified  as  follows. 
Consider  the  second-order  monic  vector  differential  equation 


q(t)  +  (Aj  +  AAi)q(t)  +  (Aj  +  AAj)q(t)  =  (B\  +  Afli) u(f)  +  D\  iv(t), 
y(t)  =  Ci  q(t)+C7q(t), 


(2a) 

(2b) 


Uncertain  linear  systems 


ft 


where  g(t)  €  Rm ,  u(t)  €  Hm ,  w(t)  €  Hm,  and  y(t)  €  Rm  are  partial  state,  input,  disturbance, 
and  output,  respectively.  The  state- variable  realization  of  the  second-order  vector  differential 
equation  in  (2)  in  a  block  companion  form  is  given  by 


x(t)  =  (A  -l-  A  A)  x(t)  +  (B  +  A  B)  u(t)  +  D  w(t),  (3a) 

y(t)  =  Cx{t),  (3b) 


where 


.  ,  _  V/  V/ 

“  -AA2  -A.4, 


C’=[C’2,C1], 


with  E  =  [-Bj  1  AA2,  -Bl  1  Aj4i],  F  =  Bll  ABi,  and  G  =  B11  D\  assuming  det(Bi)  ^  0. 
Obviously,  the  system  in  (3)  satisfies  the  matching  conditions  (i-iii)  provided  that  ||F||  <  1  or 
2I  +  F  +  F7  >0. 


Remark  1.  In  general,  if  the  uncertain  linear  system  in  (1)  satisfies  the  matching  conditions 
(i-iii),  the  matrices  E,  F,  and  G  can  be  obtained  from  AA,  A B,  and  D,  respectively,  using  a 
technique  based  on  the  singular  value  decomposition  (SVD)  [11],  I 


3.  GUARANTEED  DISTURBANCE-ATTENUATION  ROBUST-STABILIZING 
CONTROLLERS  WITH  OPTIMAL  PERFORMANCE  FOR  MATCHED  SYSTEMS 


Consider  the  following  matched  uncertain  linear  system: 

x(t)  =  (A+ BE)x(t)  +  (B+ BF)u(t)  +  BGw(t),  (4a) 

y(t)  =  Cx(t).  (4b) 


Suppose  that  the  only  information  about  the  uncertainty  matrices  in  (4)  is  that  their  matrix 
norms  are  bounded  by 

||f?||  <  a  and  ||F||</?<1.  (5) 

The  following  theorem  guarantees  that  a  disturbance-attenuation  robust-stabilizing  controller 
with  optimal  performance  exists  for  the  matched  uncertain  linear  system  in  (4)  having  the  con¬ 
straints  in  (5). 

Theorem  1.  Consider  the  matched  uncertain  linear  system  in  (4)  with  the  norm-bounded 
uncertainty  matrices  described  in  (5).  Let  6  >  0  be  any  given  disturbance-attenuation  constant 
and  Q  any  given  symmetric  positive-definite  (SFD)  matrix.  With  the  selection  of  positive  scalars 
£  i  >  0  and  £2  >  0  satisfying 


£i  < 


1-/? 

a 


and 


£  2  < 


(1  -  0  -  £i  a)i 


(G) 


(6) 


there  always  exists  a  SPD  solution  P  for  the  following  Riccati  equation: 

AtP  +  PA-  PB  [(1  — /?  — £ja)/—  ^GGrl  BT  P  +  —  /  +  —  CT  C  +  Q  =  0.  (7) 

l  0  J  £1  £2* 

Then,  a  disturbance-attenuation  robust-stabilizing  control  law  is  given  by  u(t)  =  K  x(t),  where 
K  =  -7  Br  P  with  1  >  1/2.  That  is,  the  closed-loop  system  matrix  Ae  =  A  +  B  E  +  (B  + 
B  F)  K  is  asymptotically  stable  and  the  Hoo-norm  of  the  closed-loop  transfer  functon  matrix 
H(s)  =  C(s  1  -  Ac)~l  D  (here,  D  =  BG)  is  less  than  or  equal  to  6  for  all  admissible  uncertainty 
matrices  E  and  F  in  (5).  Furthermore,  the  state-feedback  control  law  u(t )  =  —7  BT  P  x(t)  with 
7  >  1/(1  -  /?)  js  also  optimal  with  respect  to  the  following  performance  index: 


J=\J“[xT(t)Qx(t)  +  uT(t)Ru(t)\dt, 


(8a) 


ft 


Jr 


• _  •  •  •  •  • 
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where 

R=-I>  0  and  Q  =  -AT  P  -  P  A  +  P  B  ft"1  BT  P  >  0  (8b) 

7 

with  A  =  A  +  B  E  and  B  =  B  +  BF. 

Proof.  With  the  selection  of  £i  >0  and  ej  >  0  satisfying  (6),  it  is  easy  to  see  that  there  always 
exists  a  SPD  solution  P  to  the  ARE  in  (7)  [12].  To  show  the  robust  stabilization,  we  define 

Qc  =  —A'  P  —  P Ac.  (9a) 

Then 

Qc  =  -AT  P  -  PA-  E?  BT  P  -  PBE  +  yPB(2I  +  FT  +  F)BT  P.  (9b) 

From  (7),  it  follows  that 

Qc  =  PB[(2-r-l  +  0)I  +  y(FT +  F)]  BT  P  +  eiaPBBT  P 

+  —  I- EtBtP-PBE+^-PBGGtBtP  +  ct  c  +  q 

£i  6  £26 

>  <27-1)0 -nPBB'P+yZpBB-JE^yZrBE-JEi) 

+^PDDt  P+—CtC  +  Q.  (9c) 

0  £  2  v 

Hence 

Qc>e-rPDDT  P+-LctC  +  Q>0  for  ||F||<0<1  and  7  >  i,  (9d) 
0  £2  o  1 

or 

Qc>^-PDDr  P  +  -^CTC>  0  for  ||F||  <  0  <  i  and  7  >  i.  (9e) 

d  £jd  i 

Thus,  based  on  Lyapunov  stability  theory  [12],  Ae  is  asymptotically  stable  for  ||F||  <  /?  <  1  and 
7  >1/2- 

To  show  the  disturbance  attenuation,  we  utilize  the  equality  in  (9a)  and  the  inequality  in  (9e) 
as  follows: 

(-ju  l  -  Aef  P  +  P(jw  I  -  Ae)  -  %  P  D  Dt  P  -  -4cr  C  >  0  (10a) 

0  Ej  0 

for  all  w  €  Ti.  Now,  we  define  4>(jw)  =  (jw  I  —  A«)-1,  and  premultiply  DT  4>T(—jui)  and 
postmultiply  <j>(ju>)  D  to  the  inequality  in  (10a)  to  obtain 

Dt  P<t> (jw)  D  +  DT  PD-^-Dt  <t>T(-ju)  PDDt  P<t>(jw)  D 

0 

— —.Dt  <t>T{-ju)  CT  C  D>  0.  ( 10b) 

£2  « 

Then,  we  complete  a  square  term  as  follows: 

(yj/  -  DT<i,T{-ju,)PD]j  yil  -y^fDT  *TU")PDj  >  0.  (10c) 

Thus,  from  (10b)  and  (10c)  we  obtain 

7-  /  >  7^7  DT4>T(-j*)CTC<t>(ju)  D  =  — ~  HT(-ju>)  H(ju>).  (10d) 

£j  EjO  Ej  0 

Hence,  ||//(jw)||  <  6  for  all  w  €  ft- 
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To  show  the  robust  performance,  we  let  A,  B ,  R,  and  Q  be  defined  as  in  (8).  From  (9b) 
and  (9c),  we  have 

Q  =  -(A  +  BEf  P-P(A  +  BE)  +  tP(B+BF){B  +  BF)t  P 

>  PB[{I3-  1  )I  +  y{I  +  F)(I  +  F)T]BT  P+^-PDDt  P+-^CtC  +  Q. 

6  £2  0 

Since  (/  +  F)(/  +  F)T  >  (1  —  0)2  I  when  ||F||  <  0  <  1,  we  have  Q  >  0  for  7  >  1/(1  -  0).  Hence, 
the  state-feedback  control  law  u(t)  =  —7  BT  P  x(t)  with  7  >  1/(1  —  0)  is  optimal  [12]  for  the 
system  in  (4)  with  respect  to  the  quadratic  performance  index  in  (8).  I 

Remark  2.  The  Riccati  equation  in  (7)  is  constructed  to  account  for  robust  stability  and 
disturbance  attenuation  for  the  matched  uncertain  system.  If  there  is  no  system  uncertainty  (i.e., 
a  =  0  and  0  =  0)  and  disturbance  attenuation  is  not  required  (i.e.,  6  — >  00),  the  augmented  ARE 
in  (7)  reduces  to  an  ordinary  ARE  which  arises  in  the  linear  quadratic  regulator  problem  [12]. 
We  assume  Q  >  0  to  facilitate  the  proof;  however,  if  (A,C)  is  observable,  this  assumption  can 
be  relaxed  to  Q  >  0.  With  the  robust  control  law  u(f)  (=  —7  BT  P  x(t)  for  7  >  1/(1  -  0) 
and  P  >  0  being  the  solution  of  the  ARE  in  (7))  as  proposed  in  Theorem  1,  the  quadratic 
performance  index  J  in  (8),  which  is  the  compromise  of  the  weighted  state  energy  and  the 
weighted  control  energy,  can  be  minimized.  Therefore,  the  robust  control  law  u(t)  is  also  optimal 
and  provides  the  closed-loop  system  with  the  gain  margin  of  1/2  to  00  and  the  phase  margin 
of  at  least  60°  [12].  Moreover,  the  ARE  based  state-feedback  and  output-feedback  control  laws 
derived  in  [10]  provide  robust  stability  and  disturbance  attenuation  for  an  uncertain  linear  system 
with  AA  /  0  but  A B  =  0;  whereas,  our  ARE  based  state-feedback  control  law  provides  robust 
stability  and  disturbance  attenuation  for  an  uncertain  system  with  both  AA  ^  0  and  A B  /  0 
and,  also,  gives  an  additional  feature  (i.e.,  robust  performance)  for  the  same  uncertain  system. 
Furthermore,  due  to  the  simplicity  of  selecting  the  tuning  parameters  £1  and  £2  satisfying  (6),  the 
proposed  aproach  can  more  easily  determine  a  robust  control  law  for  matched  uncertain  system 
by  solving  the  ARE  in  (7)  than  the  methods  in  [4,10,13].  I 

Corollary  1.  Consider  the  matched  uncertain  linear  system  in  (4)  with  the  norm-bounded 
uncertainty  matrices  described  in  (5).  Let  6  >  0  be  any  given  disturbance-attenuation  constant 
and  Q  any  given  SPD  matrix  and  h  >  0  a  prescribed  degree  of  stability  [12].  Let  £1  and  £2  be 
any  positive  scalars  satisfying  (6),  and  P  be  the  SPD  solution  of  the  ARE: 

{A  +  hI)T  P  +  P(A  +  hi)-  PB  [(1  -0-£ia)I-e-jGGT] 

xBTP+-I+^-1CTC  +  Q  =  0.  (11) 

£1  £2  V 

Then,  a  disturbance-attenuation  robust-stabilizing  control  law  with  the  attenuation  constant  6  is 
given  by  u(t)  =  K  x(t),  where  K  =  —7  BT  P  with  7  >  1/2.  furthermore,  the  closed-loop  system 
matrix  Ae  =  A  +  B  E  +  (B  +  B  F)  K  has  a  prescribed  degree  of  stability  h  [12]  for  all  admissible 
uncertainty  matrices  E  and  F  in  (5).  I 

Now  we  consider  the  matched  uncertain  linear  system  in  (4)  with  structured  uncertainty 


matrices  E  G  72mxn  and  F  g  described  by 

k 

E  = 

»=1 

with 

|e.l  < 

and 

(12a) 

*=1 

with 

t/.  I  <  fi, 

(12b) 

respectively,  where  e,  and  f,  are  uncertain  parameters,  and  E,  and  F)  are  known  constant  matrices 
with  each  matrix  may  having  rank  greater  than  one.  Applying  the  SVD  method  [11]  to  the 
matrices  E,  and  F),  we  can  decompose  each  £,  and  F)  as 

E,  =  T,  Uj  and  F,  =  K  ttf , 

where  T,,  t/,- ,  VJ,  and  W,  are  weighted  unitary  matrices  with  appropriate  dimensions. 


(12c) 
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lb  derive  the  disturbance-attenuation  robust-stabilizing  controllers  for  the  matched  system 
in  (4)  with  the  structured  uncertainty  matrices  described  in  (12),  we  define  symmetric  positive- 
semidefinite  matrices  T  €  7Zmxm,  U  €  'R.nxn,  and  V  6  as  follows: 


e,  U,  UT,  (13a) 

i=  1  1=1 

V  =  lY.f'(V‘V?  +  W*W?'>'  (13b> 

1  .= l 

with  the  matrices  T),  Ui,  Vi,  and  Wj-  as  in  (12).  It  can  be  shown  that  2V  +  F  +  FT  >  0.  Also, 
from  the  matching  condition  (ii),  we  require  2/  +  F  +  F7  >  0.  As  a  result,  we  assume  that 

I-V>  0.  (13c) 

The  following  theorem  guarantees  that  a  disturbance-attenuation  robust-stabilizing  controller 

with  optimal  performance  exists  for  the  matched  uncertain  linear  system  in  (4)  with  the  structured 
uncertainty  matrices  in  (12). 

Theorem  2.  Consider  the  matched  uncertain  linear  system  in  (4)  with  the  structured  uncertainty 
matrices  described  by  (12).  Let  6  >  0  be  any  given  disturbance-attenuation  constant  and  Q  any 
given  SPD  matrix.  With  the  selection  of  positive  scalars  C\  >  0  and  ti  >  0  satisfying 

r  ^  I- <J,max(l^)  ,  r  [I- 0max(lO  —  fl  ‘’mw  (T)]6  .  . 

£i  <  — TrT  411(1  £*  - Is — Tr \ - >  l14) 

(T)  ffm»x  (C) 

there  always  exists  a  SPD  solution  P  for  the  following  Riccati  equation: 

At  P+PA-PB(l-V-e1T-~GGT)  BtP  +  -U  +  -^CtC  +  Q  =  0,  (15) 

'  v  *  €\  £2  0 

where  the  matrices  T,  U,  and  V  are  defined  in  (13).  Then,  a  disturbance-attenuation  robust- 
stabilizing  control  law  with  the  attenuation  constant  6  is  given  by  u(t)  =  K  x(t),  where 
K  =  —7  Bt  P  with  7  >  1/2.  Furthermore,  the  state-feedback  control  law 


u(<)  =  —7  Bt  P  x(t)  with  7  >  - - - 

1  —  <rm 


'  '  '  '  -  l-<rm„(K) 

is  also  optimal  with  respect  to  the  quadratic  performance  index  as  defined  in  (8). 

Proof.  With  £1  >  0  and  e2  >  0  satisfying  (141,  it  is  easy  to  see  that  there  always  exists  a  SPD 
solution  P  to  the  ARE  in  (15)  [12].  Define  Qc  as  in  (9a).  From  (15),  it  follows  that 

Qc  =  PB  [(27  -  1)  /  +  V  +  7  (FT  +  F)]  Bt  P  a  t j  P  B  T  BT  P  +  U 
-ETBTP-PBE  +  ^PBGGTBTP+—CTC  +  Q. 

0  1 2  v 


2 v  +  ft  +  f  =  [fi(Vi  V?  +  w,  W?)  +  MVi  w?  +  Wi  v?)} 


>  ^2  l/<l  (K  ±  W.XK  ±  Wif  >  0 


f 
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eiPBTbTP  +  ~U~ETBTP-PBE 
£1 


=  Xj[e.'  (e\P  BTif?  BT  P  +  -ei(UiTT  BT  P  +  PBTiU?)  J 

>  J2 m  (y^p BTi±~jrx  Ui)  (v^F BTi±^)ff u,s)  -  °- 


It  follows  that 


Q c  >  PB[(2-,-l)I  +  V -2-rV]BT  P+£-lpBGGT  BT  P+—  CTC  +  Q 

6  e2  6 

=  (2-r-\)PB(I-V)BT  P  +  £-lpBGGT  BT  P+-L-CtC+Q. 

6  £jS  ^ 


Hence, 


Qc  >  t  P  D  Dt  P  +  — p  CT  C  +  Q  >  0  for  I  —  V  >  0  and  7  >  i 
d  e2d  ~  2 

Thus,  based  on  Lyapunov  stability  theory  [12],  Ae  is  asymptotically  stable  for  /  -  V  >  0  and 

7  >1/2. 

The  proofs  for  disturbance  attenuation  and  the  optimality  condition  are  similar  to  those  in 
Theorem  1  and  hence  omitted.  g 

Remark  3.  Note  that  the  robust  control  law  obtained  in  Theorem  1  is  more  conservative  than 
that  obtained  in  Theorem  2  due  to  different  uncertainty  structures.  In  general,  the  control  gain 
obtained  in  Theorem  1  is  larger  than  that  obtained  in  Theorem  2.  | 

4.  DISTURBANCE-ATTENUATION  ROBUST-STABILIZING  CONTROLLERS  FOR 

MISMATCHED  SYSTEMS 


•  • 


Consider  the  following  mismatched  uncertain  linear  system  described  by 

*(t)  =  (A  +  AA)  i(f)  +  (B  +  AB)ti(t)  +  Dw(t),  (16a) 

y(<)  =  C  x(t).  (16b) 

Suppose  that  the  only  information  about  the  uncertainty  matrices  AA  6  Unxn  and  AH  6  Unxm 
in  (16)  is  that  the  matrix  norms  are  bounded  by 

||AA||<a  and  ||AH||  <  0.  (17) 

The  following  theorem  will  be  utilized  to  find  a  disturbance-attenuation  robust-stabilizing  con¬ 

troller  for  the  mismatched  uncertain  system  in  (16)  with  the  constraints  in  (17). 

Theorem  3.  Consider  the  mismatched  uncertain  system  in  (16)  with  the  norm-bounded  uncer¬ 
tainty  matrices  described  in  (17).  Let  6  >  0  be  a  given  disturbance-attenuation  constant  and  Q 
a  given  SPD  matrix.  Suppose  that  there  exist  positive  scalars  C\  >  0,  Ci  €  (0, 2/0) ,  and  s3  >  0 
such  that  the  Riccati  equation 

atp  +  pa-p^i-£jI^  BflT-(Eia  +  l)  1~jddT 

xP+j;I+7r6cTc+Qzz0  (18) 


t 
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has  a  SPD  solution  P.  Then,  a  disturbance-attenuation  robust-stabilizing  control  law  is  given  by 
u(t)  =  K  *(<),  where  K  =  -7  BT  P  with  7  satisfying  either 

»  jarajr  !>»  <19> 

That  is,  the  closed-loop  system  matrix  Ac  =  A  +  A  A  +  (B  +  A  B)  K  is  asymptotically  stable  and 
the  Hoo -norm  of  the  closed-loop  transfer  functon  matrix  H(s)  =  C  (si  —  j4e)_1  D  is  less  than  or 
equal  to  S  for  all  Omissible  uncertainty  matrices  A  A  and  A B  in  (17). 

Proof.  Suppose  that  the  Riccati  equation  in  (18)  has  a  SPD  solution  P.  Define  Qc  as  in 
Theorem  1.  From  (18),  it  follows  that 

Qc  =  P  [^27-1  +  ^)  BBT  +  ^-I  +  7BABt  +  7*BBt  P 

+  (e1aPP+  —  I-AATP-PAA')+^PDDTP+-^CTC  +  Q. 

V  £1  )  °  £ao 


l72fi0BBT  +  —  I  +  7BABt  +  7ABBt 
2^2 

>(71/^b+^ab)(,^b+-?5L5ab)t>o 


£iC,PP  +  ~I-AATP-PAA 


we  obtain  the  following  inequality: 

Qc  >  (27- l  +  ^-  272£2^  PBBt  P  +  ^-PDDt  P+-^CtC+Q 
=  (27  —  1)  l-£l^(27+l)  pbbt  P+^-PDDt  P  +  —-CtC  +  Q. 

[2  0  €3  0 

If  7  satisfies  either  inequality  in  (19),  which  implies 

(27-1)  l-^(27+l)  >0,  then  Qc>^PDDTP^^-CTC  +  Q>  0. 

2.  0  £30 


Thus,  based  on  Lyapunov  stability  theory  [12],  the  obtained  controller  u(t)  stabilizes  the  mis¬ 
matched  system  in  (16)  with  the  constraints  in  (17). 

The  proof  for  ||//||oo  <  6  is  similar  to  that  in  Theorem  1  and  hence  omitted.  I 

Remark  4.  The  parameter  £3  in  (18)  is  restricted  to  be  in  the  range  of  (0,2//?)  such  that  the 
term  (1  —  £3  0/2)  in  (18)  is  greater  than  zero.  I 

Now  we  consider  the  uncertain  linear  system  in  (16)  with  structured  uncertainty  matrices 
Aj4  €  Kn*n  and  AB  €  Kn*m  described  by 

* 

Ai4  =  ^a,  j4j  with  |a^  |  <  dj,  and  (20a) 

i=i 

1 

AB  =  £bjBi  with  |M<6j,  (20b) 


Uncertain  linear  eyitema 


75 


respectively,  where  a,-  and  6,  are  uncertain  parameters,  and  Ai  and  Bi  are  known  constant 
matrices  with  each  matrix  may  having  rank  greater  than  one.  Applying  the  SVD  method  [11]  to 
Ai  and  we  can  decompose  each  A,  and  Bi  as 

A,  =  Ti  Uj  and  B,  =  Vi  Wf,  (20c) 

where  Ti,  £/*,  Vi,  and  Wi  are  weighted  unitary  matrices  with  appropriate  dimensions. 

To  derive  the  disturbance-attenuation  robust-stabilizing  controllers  for  the  system  in  (16)  with 
the  structured  uncertainty  matrices  described  by  (20),  we  define  symmetric  positive-semidefinite 
matrices  Tg  Knxn ,  U  €  Hnxn ,  V  eHnxn,  and  W  €  Kmxm  as  follows: 

T  =  £  a,  7)  7f,  UlJ^aiUil/?,  (21a) 

*=1  1=1 

v  =  l  Tbivi  v;T,  w  =  \  T  h  Wi  wj,  (21b) 

i  t=i  i=i 

with  the  matrices  7),  Ui,  Vi,  and  Wi  as  in  (20).  The  following  theorem  will  be  utilized  to  find  a 
disturbance-attenuation  robust-stabilizing  controller  for  the  mismatched  uncertain  system  in  (16) 
having  the  constraints  in  (20). 

Theorem  4.  Consider  the  mismatched  uncertain  linear  system  in  (16)  with  the  structured 
uncertainty  matrices  described  in  (20).  Let  S  >  0  be  a  given  disturbance-attenuation  constant 
and  Q  a  given  SPD  matrix.  Suppose  that  there  exist  positive  scalars 


£l  >0,  52  G 


such  that  the  Riccati  equation 


(°’<WW))’ 


and  £3  >  0 


At  P  +  P  A-  P  (b  Bt  -  £l  T  -  £iBW  Bt  -  j-V  -  £-j  D  DT^ 

x  P  +  —  +  -^7  CT  C  +  Q  =  0  (22) 

5l  53  0 

has  a  SPD  solution  P,  where  T,  U,  V,  and  W  are  defined  in  (21 ).  Then,  a  disturbance-attenuation 
robust-stabilizing  control  law  with  the  attenuation  constant  6  is  given  by  u(t)  =  K  x(t),  where 
K  —  —7  BT  P  with  7  satisfying  either 

Proof.  Suppose  that  the  Riccati  equation  in  (22)  has  a  SPD  solution  P.  Define  Qe  as  in 
Theorem  1.  From  (22),  it  follows  that 

Qc  =  P\(2-r-l)BBT  +£2BWBt  +  —  V  +  yB£sBT  +i&BBt  P 

L  52 

+  (elPTP+—U-AATP-PAA)+^-PDDTP+-^-1CfrC  +  Q. 

\  5i  53  b 


4j7£3BWBt+  —  V  +  7BA£r+7ABBr 
52 


>X>I  (yy/K1BWi±-^V?j(iV2r3BWi±-±=V^  >0 
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and 


e1PTP+—U-&AI  P-PA.A 
£1 


k  T 

•  £  w  ±  ^  u>)  (y^pTi  ±  >o, 


we  obtain  the  following  inequality: 

<2C  >  PB[(27-l)/  +  f2^-  472£2W^]5TP  +  ^PDDTP+-i-CTC  +  Q 

0  £36 

=  (27-l)PB[/-e2(27  +  l)W]  flT  P  +  ^  P D DT  P  +  CT  C  +  Q. 

0  £30 


If  7  satisfies  either  inequality  in  (23),  which  implies 

(27  -  1)  [/  -  £2  (27  +  1)  W]  >  0,  then  Qc  >  ^  P  D  DT  P  +  ~  CT  C  +  Q  >  0. 

0  £36 

Thus,  based  on  Lyapunov  stability  theory  [12],  the  obtained  controller  u(t)  stabilizes  the  mis¬ 
matched  system  in  (16)  with  the  constraints  in  (20). 

The  proof  for  ||P||oo  <  6  is  similar  to  that  in  Theorem  1  and  hence  omitted.  | 

Remark  5.  The  introduction  of  tuning  parameters,  £1,  £2,  and  £3  in  (18)  and  (22),  makes  the  pro¬ 
posed  approach  more  flexible  in  obtaining  disturbance-attenuation  robust-stabilizing  controllers. 
For  instance,  consider  the  following  Riccati  equation: 

AtP+PA-P  (bBt-^DDt'J  P  +  CtC  =  0,  (24) 

which  is  the  ARE  for  the  standard  control  problem  (i.e.  the  control  effort  u(t)  is  included 
in  the  controlled  output  y(t))  in  [7].  Now,  if  there  exists  a  P  >  0  satisfying  (24)  with  ( A,C ) 
observable,  then  u(t)  =  —  (1/2)  BT  P  x(t)  can  be  interpreted  as  a  disturbance-attenuation  con¬ 
troller  for  the  Hoo  control  problem  associated  with  (16)  (i.e.  u(f)  is  not  included  in  y(<)).  It  is 
seen  tha,  (24)  corresponds  to  a  special  case  of  (18)  or  (22)  (when  A  A  =  0  and  A  B  =  0)  with 
£3  =  1/6  and  Q  =  0.  Hence,  by  adjusting  the  tuning  parameter  £3,  the  possibility  of  finding  a 
SPD  solution  for  (22)  is  greatly  enhanced  over  that  for  (24).  Also,  it  should  be  noted  that  the 
inequality  in  (23)  gives  an  explicit  bound  for  which  the  control  gain  is  allowed  to  vary  without 
affecting  robust  stability  and  disturbance  attenuation  of  the  closed-loop  system.  | 


5.  ILLUSTRATIVE  EXAMPLES 


* 


•  • 


Example  1.  Consider  a  version  of  the  pitch-axis  model  for  the  AFTI/F-16  flying  at  3000  ft.  and 
Mach  0.6  [4,13,14].  The  equations  of  motion  are  represented  in  the  state-space  form  as 


*(f)  =  (A  +  AA)x(t)  +  {B  +  AB)u(t)  +  Dw(t), 
y(t)  =  C  x(t), 

where  the  nominal  system  are  described  by 


'  0 

1 

0  ‘ 

0 

0  ‘ 

A  = 

0 

-0.87 

43.22 

,  B  = 

-17.25 

-1.58 

0 

0.99 

-1.34 

-0.17 

-0.25 

O 

II 

'  1  0 

1  0 

1  ' 
0 

J 

•  • 


•  •  • 


•  • 
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and  the  structured  uncertainty  matrices  are  described  by 


‘  0 

0 

0  ' 

'  0 

0  ' 

AA  = 

0 

a2 

A  B  = 

k 

h 

0 

03 

a*  J 

bz 

k 

with  Iml  <  0.7,  |n2|  <  35,  |a3|  <  0.7,  |a4|  <  1.05,  |6i|  <  2,  |62|  <  0.2,  |63|  <  0.02,  and  |64|  <  0.03. 

Note  that  this  system  is  matched  and  the  structured  uncertainty  matrices  can  be  expressed 
as  AA  =  B  E  and  A B  =  B  F ,  where 


[0  -0.0618  cu  +  0.3907  a3  -0.0618  a2  +  0.3907  a4  1 

h  ~  0  0.0420  ai  -  4.2657  a3  0.0420  a2  -  4.2657  a4  J 

and 

_  -0.0618  bi  +  0.3907  i3  -0.0618  f>2  4-  0.3907  fc4  ' 

*  ~  0.0420  6i  -  4.2657  63  0.0420  b2  -  4.2657  64  1 

and  the  disturbance  matrix  can  be  written  a s  D  —  BG  with 


-0.0618  0.3907 

0.0420  -4.2657 


The  eigenvalues  of  A  Eire  —7.65,  0,  5.44  and  the  nominal  system  is  unstable.  To  find  a  disturbance- 
attenuation  robust-stabilizing  control  law  for  this  matched  uncertain  system,  we  determine  T,  U , 
and  V  as  in  (13)  and  obtain 


1.8874  -1.9219 

-1.9219  8.2777 


U  =  diag  [0,3.0508,7.1143], 


and 


0.17472  -0.04797 
-0.04797  0.20393 


Set  the  disturbance-attenuation  constEint  6  =  1  and  choose  Q  =  /,  e\  =  0.04  €  (0,0.086),  and 
e2  =  0.01  €  (0,0.022).  The  Riccati  equation  in  (15)  has  a  SPD  solution 


P  = 


122.72  0.8920  3.1551 

0.8920  0.5816  -0.0804 

3.1551  -0.0804  54.211 


Then,  from  Theorem  2,  a  disturbance-attenuation  robust-stabilizing  control  law  with  6  =  1  can 
be  constructed  as  u(t)  =  K  *(<),  where 


K  =  -7  Bt  P  =  7 


15.924  10.019  7.8291  ' 
2.1982  0.8988  13.426 


with  7  >  -. 


Furthermore,  the  state-feedback  control  law 


u(t)  = -f  BT  P x(t)  with 


1 

7-  i-<w(n 


1.3149 


is  optimal  with  respect  to  the  quadratic  performance  index  in  (8). 

To  guarantee  that  the  closed-loop  system  has  a  prescribed  degree  of  stability  h  =  1,  we  set  6, 
Q,  £x,  e2  as  before  and  replace  A  by  A  4- 1  to  solve  the  ARE  in  (15)  for  P.  Then,  a  disturbance- 
attenuation  robust-stabilizing  control  law  with  6  =  1,  which  also  guarantees  that  the  states  decay 
no  slower  than  e_l,  can  be  constructed  as  u(<)  =  K  x(t),  where 


K  =  -7  Bt  P  =  7 


33.018  10.236  4.7750 
-6  4293  0.8007  20.907 


with 


7  > 


1^ 

2 


•  •••••• 
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When  the  requirement  of  disturbance  attenuation  is  relaxed,  i.e.  S  — *  oo,  a  robust-stabilizing 
control  law  u(t )  =  K  x(t)  =  —y  BT  P x(t)  for  the  matched  system  is  determined  by  solving  the 
ARE  in  (15)  for  P  with  Q  —  l  and  £\  =  0.04  as  before.  The  feedback  gain  is  given  by 


K  =  -yBT  P-y 


5.6870  6.6475 
-0.1324  0.7230 


10.092 

3.2596 


with 


*4 


Note  that  even  with  A B  ^  0,  the  obtained  control  gain  is  smaller  in  magnitude  than  those 
obtained  in  [4,13]  for  the  same  uncertain  system  but  with  A B  =  0.  Moreover,  the  proposed 
method  is  easier  to  use  in  obtaining  a  robust-stabilizing  control  law  than  those  in  [4,13],  beacause 
only  one  Riccati  equation  needs  to  be  solved  for  the  proposed  approach.  I 

Example  2.  The  dynamics  of  a  helicopter  in  a  vertical  plane  for  an  airspeed  range  of  60-170 
knots  are  given  in  [4,15],  There  are  four  state  variables — xi  =  horizontal  velocity  (knot/sec), 
x2  =  vertical  velocity  (knot/sec),  X3  =  pitch  rate  (deg/sec),  and  x4  =  pitch  angle  (deg) — and 
two  control  variables — tij  =  collective  pitch  control  and  u2  =  longitudinal  cyclic  pitch  control. 
In  the  airspeed  range  of  60  knots  to  170  knots,  significant  changes  occur  only  in  element  032,  <134, 
and  621  •  For  this  range  of  operating  conditions, 


A  = 


'  -0.0366 

0.0271 

0.0188 

-0.4555  1 

0.4422 

0.1761 

0.0482 

-1.01 

0.0024 

-4.0208 

3.0447 

-7.5922 

0.1002 

0.2855 

-0.707 

1.3229 

-5.52 

4.99 

0 

0 

1 

0 

0 

0 

D=  [0,0,0,  l]7 


C=  [0,1, 0,0], 


'  0 

0 

0 

0 

0 

0  ' 

0 

0 

0 

0 

,  A  B  = 

S21 

0 

0 

t*32 

0 

T34 

0 

0 

0 

0 

0 

0 

L  0 

0 

AA  = 


with  |r32|  <  0.2192,  |r34|  <  1.2031,  and  |*2I|  <  2.0673.  Define  T,  U,  V,  and  W  as  in  (21)  and 
obtain 


7  =  diag  [0,0, 1.4223,0], 
V'  =  diag  [1.03365,0], 


U  =  diag  [0, 0.2192,0,1.2031], 
W  =  diag  [0,1.03365,0,0]. 


Set  the  disturbance-attenuation  constant  6  =  0.5  and  choose  Q  =  /,  £1  =  1,  £2  =  0.25  and 
£ 3  =  0.25,  the  Riccati  equation  in  (22)  has  a  SPD  solution 


P  = 


9.9891  -0.6427  -1.2810  -11.2650 

-0.6427  1.0287  0.8892  2.0922 

-1.2810  0.8892  1.2521  3.4268 

-11.2650  2.0922  3.4268  19.4367 


Then,  from  Theorem  4,  a  disturbance-attenuation  robust-stabilizing  controller  can  be  constructed 
as  u(t)  =  K  x(t)  =  -7  Bt  Px(t),  where 


K--yBrP- 


-9.5318  2.0603  4.7707  17.5269 
-0.2459  3.4864  0.7284  0.7682 


^L__-i  =  l.2093  >7>i. 


with 
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To  show  the  flexibility  of  the  proposed  method  due  to  the  introduction  of  the  tuning  pa¬ 
rameters,  we  let  AA  =  0  and  A B  -  0  (i.e.,  T  =  0,  U  =  0,  V  =  0,  and  W  =  0),  and  set  the 
disturbance-attenuation  constant  6  =  0.1.  The  ARE  in  (24)  which  is  now  identical  to  (22)  with 
Q  =  0  and  £3  =  1/6  =  10  does  not  have  a  SPD  solution;  however,  with  Q  =  0  and  by  adjusting 
£3  =  0.25,  the  ARE  in  (22)  has  a  SPD  solution.  Hence,  the  desired  disturbance-attenuation 
state-feedback  control  gain  with  6  =  0.1  is  given  by 


A' =7 


-0.0033  -2.1201  0.2444  0.4382 

0.0063  5.8232  0.0755  -0.3804 


for 


Thus,  the  introduction  of  the  tuning  parameters  indeed  enhaces  the  flexibility  of  the  proposed 
method  in  finding  the  disturbance-attenuation  robust-stabilizing  controllers.  Note  that  the  above 
comparison  does  not  imply  that  the  solution  in  (24)  is  conservative  because  (24)  is  originally 
derived  for  the  standard  control  problem  with  u(t)  inc’uded  in  the  controlled  output  y(t). 
However,  when  dealing  with  the  disturbance  attenuation  control  problem  in  (16)  with  u(f)  not 
included  in  y(t),  (22)  does  lead  to  better  disturbance  attenuation  (smaller  6)  than  (24)  due  to 
the  introduction  of  the  tuning  parameters  in  (22).  I 

Remark  6.  While  the  introduction  of  tuning  parameters  provides  additional  flexibility,  the  ap¬ 
plication  of  Theorem  4  to  a  given  mismatched  uncertain  linear  system,  in  general,  may  not  always 
lead  to  a  robust  control  law.  However,  in  our  other  simulation  examples,  we  have  successfully 
determined  various  robust  control  laws  via  appropriate  adjustment  (i.e.,  successive  reduction)  of 
the  tuning  parameters,  £t,  £2,  and  £3  in  (22),  without  numerical  problems.  | 


6.  CONCLUSION 


Based  on  the  LQR  theory  and  Lyapunov  stability  theory,  new  disturbance-attenuation  robust- 
stabilizing  controllers  have  been  developed  for  matched  and/or  mismatched  uncertain  linear  sys¬ 
tems.  It  has  been  shown  that  dynamic  systems,  described  by  second-order  vector  differential  equa¬ 
tions,  often  satisfy  the  matching  conditions  and  that  disturbance-attenuation  robust-stabilizing 
controllers  (with  optimal  performance  if  ||Afl||  <  1/2)  always  exist  for  matched  uncertain  linear 
systems  which  contain  structured  or  norm-bounded  uncertainty  matrices.  For  mismatched  un¬ 
certain  linear  systems  ,  two  theorems  have  been  developed  for  finding  disturbance-attenuation 
robust-stabilizing  controllers.  These  disturbance-attenuation  robust-stabilizing  control  laws  can 
be  easily  constructed  from  the  symmetric  positive-definite  solution  of  the  augmented  Riccati  equa¬ 
tion.  Also,  the  proposed  approach  is  more  flexible  than  some  existing  methods  in  the  sense  that 
additional  tuning  parameters  (such  as  e,  7,  and  h,  etc.)  have  been  introduced  in  the  derivations 
to  achieve  robust  stabilization,  robust  performance,  and  disturbance  attenuation  for  uncertain 
linear  systems.  Two  practical  examples  have  been  presented  to  illustrate  the  results. 
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